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We obtain all "regularization parameters" (RP) needed for calculating the gravitational and 
electromagnetic self forces for an arbitrary geodesic orbit around a Schwarzschild black hole. These 
RP values are required for implementing the previously introduced mode-sum method, which allows 
a practical calculation of the self force by summing over contributions from individual multipole 
modes of the particle's field. In the gravitational case, we provide here full details of the analytic 
method and results briefly reported in a recent Letter [Phys. Rev. Lett. 88, 091101 (2002)]. In the 
• electromagnetic case, the RP are obtained here for the first time. 
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I. INTRODUCTION 



This is the second in a series of papers aimed to establish a practical calculation scheme for the self force acting 
on a point particle in orbit around a black hole. This scheme — referred to as the "mode-sum method" — stems from 
^ the general regularization prescription by Mino, Sasaki, and Tanaka (MST) and Quinn and Wald (QW) ||]. In 
CN| ' effect, the mode-sum method re-formulates the MST-QW general result in the language of multipole modes, thereby 
making it accessible to standard numerical treatment. In practice, the application of this method involves two basic 
parts: (i) calculation of the "full" modes of the force, through numerical integration of the decoupled field equations; 
, and (ii) analytical derivation of certain parameters (whose values depend on the orbit under consideration) called the 
• "regularization parameters" (RP). Previously, the explicit values of the RP were derived analytically in a few special 
\ cases of orbits in Schwarzschild spacetime — specifically, for circular and radial orbits in the scalar field case |^ , and 
for radial trajectories in the gravitational case |^,||. In these works, the RP values where calculated through a rather 
cumbersome local expansion of the (/-multipole) Green's function. The application of this technique to more general 



I orbits appears a challenging task, 
^ijpj' In a recent Letter the joint groups of Barack and Ori (BO) and Mino, Nakano, and Sasaki (MNS) devised an 
. . , alternative, more direct method for obtaining the RP. The new method is based on a multipole decomposition of the 
^ • explicit "direct" part of the force (see below). Using this method, BO and MNS were able to calculate the explicit 
k>( \ RP values for both the scalar and gravitational self forces, for any geodesic orbit in Schwarzschild spacetime. In a 
j_] ' preceding paper |^ (hereafter referred to as "paper I") BO described the full details of the new calculation technique, 
■ as applied to the toy-model of the scalar self force acting on a scalar charge. The current paper deals with the more 
interesting case of the gravitational self force on a mass particle, and provides full details of the RP derivation in this 
case. In addition, we shall derive here the RP values for the electromagnetic self force acting on an electrically charged 
particle orbiting a Schwarzschild black hole. (BO and MNS applied two slightly different methods in obtaining the 
RP; MNS describe their calculation in [||.) 

The analysis presented in this paper relies greatly on the technique and results of Paper I, to which we shall 
frequently refer the reader. Though the basic idea of our calculation is the same as in the scalar toy model of Paper 
I, some unavoidable technical complexities arise when coming to consider the gravitational or electromagnetic cases. 
In particular, one then has to consider an extension of the particle's four-velocity (which takes part in constructing 
the "direct force" — see below) off the worldline, and address the question of the RP dependence on the (non-unique) 
choice of such extension. Another, more fundamental issue, is the gauge dependence of the gravitational self force 
and its implication to the mode sum scheme (see Ref. 1^). 

Most of our manuscript will concern with the (most interesting) gravitational case. Our analytical technique is easily 
applicable to the (simpler) electromagnetic case, which we shall later consider in a separate section. The structure 
of this paper is as follows: In the rest of this introductory section we set up the physical scenario of a pointlike 
mass particle orbiting a Schwarzschild black hole, introduce MST-QW's prescription for calculating the gravitational 
self force on this particle, and review the basics of the mode-sum method. In Sec. |l| we present MST's expression 
for the direct part of the gravitational self force, and analytically process this expression to extract the information 



1 



relevant for deriving the RP. Section II] contains the heart of our calculation, namely, the derivation of all RP for 
any geodesic orbit, in the gravitational case. Section IV deals with the electromagnetic case. In Sec. we summarize 
the prescription for calculating the gravitational and electromagnetic self forces via the mode-sum method, and give 
some concluding remarks. 

Throughout this paper we use geometrized units (with G = c = 1), metric signature — h++, and the standard 
Schwarzschild coordinates t, r, 6, (p. 



A. Pointlike particle model 



We consider a pointlike particle of mass /i, moving freely in the vacuum exterior of a Schwarzschild black hole with 
mass M 3> /i. (QW Q discuss the extent to which the concept of a pointlike particle makes sense in the context of the 
radiation reaction problem.) In the limit /j, — > 0, the particle traces a geodesic z^(r) of the Schwarzschild background. 
Due to angular momentum conservation, the geodesic orbit (as well as the orbit under self-force effect) is confined to 
a plain, which, without loss of generality, we shall take as the equatorial plain, 6 = 7r/2. 

When the mass /i is finite, the particle no longer moves on a geodesic. In this case, it is useful to view the particle 
as being subject to a self force induced by its own gravitational field (treated as a perturbation over the background 
geometry). a The particle's equation of motion thus takes the form 

^xu^.^puf' = (1) 

where u" = dz'^/dr, a semicolon denotes covariant differentiation with respect to the background geometry, and 
oc 0{fi^) describes the leading-order self-force effect. In the following we shall consider the self force acting on 
the particle at an arbitrary point along its worldline, denoted by z = zq = (tg, tq, ^Oj fo) (in our setup Oq = 7r/2). We 
shall use the notation x = {t, r, 9, ip) to represent a point in the neighborhood of zo- 

We will denote the metric of the perturbed spacetime as gap + haf3, where gap is the (Schwarzschild) background 
metric and hap{o^ fi) is the metric perturbation induced by the particle. Following MST-QW, we consider the metric 
perturbation haf3 specifically in the harmonic gauge. We shall denote by haf3 the trace-reversed perturbation: 

hal3 = hap - -zgaf3- (2) 



B. Gravitational self force according to MST-QW 



MST and QW found that the gravitational self force on a particle freely moving in a vacuum spacetime can be 
formally constructed as ||l^ 

F,Zu - lim F,l,,{x), (3) 

where Jj^jj is the "tail" force, associated with the mere effect of waves scattered inside (rather than propagating 
along) the particle's past light cone. The tail force can be derived from the "tail" part of the metric perturbation, as 
defined by MST through 

F,l^{x)=fik-r'^\x)hlf^,{x). (4) 
Here, the tensor k"^^^{x) is any (sufficiently regular) extension of the tensor 



^One might attempt an alternative point of view, which regards the particle as freely moving in a "perturbed spacetime". 
This picture, though, is somewhat problematic, as the perturbed spa cetime is singular at the point-particle limit, and hence is 
not strictly defined at the particle's location. See, however, Ref. [[lO|. 
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defined at a; = zq, where m" and g"^ refer to the values of the four- velocity and the metric tensor at zq. [Later we 
shall employ a specific extension of k^^^^; note that the choice of extension does not affect the physical self force F^^^, 
though, obviously, it does affect the field F^■^{x) off the worldhne.] 

The (singular) difference between the "full" perturbation haplx) and the tail part h^^p{x) is associated with the 
instantaneous effect of waves propagating directly along the particle's light cone. This part is referred to as the 
"direct" perturbation: 

h%{x)^Kp{x)-h':^{x). (6) 
Correspondingly, we define the "direct" force as 

F2,,{x)^^ik'^P'^'{x)hf^.^,{x). (7) 

Defining also the "full" force, 

FU{x)^^xk'^P^\x)hp,,s{x), (8) 

we then have 

FZ^^ix) = Ft^,,{x) ~ F^x). (9) 

The explicit form of the direct perturbation has been derived by MST (see also |pl],^). It is given below in Eq. 
(|l3|), and serves as the starting point for our analysis. 

Note that both the direct force and the full force, which were defined above as vector fields in the neighborhood 
of Zq, diverge as x — > zq. However, their difference, yielding the tail force, admits a perfectly regular limit x — > zq, 
which, according to MST-QW, represents the physical self force. In this respect, notice also the freedom one has in 
choosing the extension k^^^'^ix), as long as this extension is regular enough and reduces to k^^'^^ in the limit a; — > zq. 
One has to make sure, though, that the same extension k°'^^^{x) is applied to both the direct and the full forces. 



C. Mode-sum method 



The mode-sum method was previously introduced as a practical method for calculating the MST-QW self 
force given in Eq. (^) ||l^ , |l3|| . The method is reviewed in Paper I; here we merely describe the basic prescription (as 
applied to the gravitational case) and introduce the relevant notation. 

In the mode-sum scheme, one first formally expands the gravitational tail force, as well as the full and direct forces, 
into multipole /-modes, in the form 



^tail,full,dir(2^) 



tail. full. dir 



(x). 



(10) 



1=0 



Here, precisely as in the scalar case, the modes -Fjaii' ^fu'u' ^^'^ ^dil obtained by decomposing (each of the vectorial 
components of) the corresponding quantities F^^^, ^f^n? ^-iid F^^^. into standard scalar spherical harmonics, and then, 
for any given multipole number I, summing over all azimuthal numbers m. It is important to emphasize here that 
the various Z-modes introduced in Eq. (10) are defined in our scheme through a scalar harmonic decomposition. In 
this regard, recall that the (full) metric perturbation in Schwarzschild spacetime is usually decomposed into tensor 
harmonic modes in actual calculations. The construction of the full-force scalar-harmonic modes ff" ^, from the full 
perturbation tensor-harmonic modes can be prescribed in a straightforward manner (as, e.g., in ||5|,|l4|). 

The basic prescription for constructing the gravitational self force via the mode-sum scheme is given by MM 



lim Ff"J^{x) - A'^L - - C"/L 

x^zo 



(11) 



where L = I + 1/2 and the (/-independent) coefficients A", B", C", and 13" are the regularization parameters (RP). 
The RP A", B", and C" may be defined by the demand that the sum in Eq. (11) would converge. Equivalently (and 
more practically) , one may define these parameters by requiring convergence of the sum 
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lim F<j^'.(a;) - A" L - B" - C"^ / L 



(12) 
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This sum then defines the fourth parameter, D". From the above definitions it is clear that the RP values may be 
derived through analysis of the direct force modes F^i^{x). 

Eq. (|ll|) constitutes a practical prescription for constructing the gravitational self force, given (i) the values of all 
necessary RP, and (ii) the full- force modes F^^^. In this paper we derive all RP for any (equatorial) geodesic orbit in 
Schwarzschild spacetime, hence setting an analytical basis for calculations of the gravitational self force for all such 
orbits. 



II. ANALYZING THE DIRECT GRAVITATIONAL FORCE 

A. Direct part of the metric perturbation 

The direct part of the trace-reversed metric perturbation was obtained by MST — see Eq. (2.27) of Ref. [Q. In the 
Appendix we process the expression obtained by MST, and bring it to the form 

hf:;{x) = 4/ie-i uliil^ + e-^pf^{x, zq). (13) 

Here, e is the spatial geodesic distance from the point x to the geodesic z{t) (i.e., the length of the short geodesic 
connecting x to the worldline and normal to it), zi denotes the intersection of this short normal geodesic with the 
worldline, and v}^ is the four-velocity parallelly-propagated (PP) from zi to x. (See Fig. [y for an illustration of the 

geometric setup described here.) The function P^^ is a regular function of x, of order 5x^ (and higher orders), where 
5x^ = x^ — Zq. The explicit form of Pi^' will not be needed in our analysis. 




\ 

FIG. 1. Geometric setup and notation: The self force is calculated at the point ^ = zq on the geodesic z{t). x is an off-worldline point 
in the neighborhood of zq, and 5x = x denotes the length of the short geodesic section connecting x to the worldline and 

normal to it, and z\ is the point where this short geodesic intersects the worldline. u°' and ii^* denote the four- velocities at zq and zi, 
respectively. u°' and (the former not shown in the sketch) are vectors at x, generated by parallelly-propagating the vectors n° and 
along the short geodesic section from zq or zi, respectively, to x. 

For later convenience, we first re-express /i^'^ in terms of the four-velocity PP from zq to x (rather than from zi to 
x), which we denote by Ua (or tt"). Both Ua and are regular functions of a;, and the difference between them is 

proportional to 5x^ (and to the Riemann tensor). Absorbing this difference in the function P'jj^^ , we may re- write the 
direct metric perturbation as 

hf^{x) = ^iiS-^'HpiL^ + S-^'^pf^{x, zo), (14) 
where S = e^, and the new function Pg^-* has the same features as -Pg^\ namely, it is a regular function, of order 5x^ . 

B. Extending the tensor k^^^^ off the worldline 

Given the above expression for the direct perturbation, the direct force is constructed as a vector field through 
Eq. (|^). In this equation, recall, k"^^^{x) is an extension off the point zq of the tensor k^^^^ defined at zq in terms 
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of u°'{zo) and g°'^{zo) [see Eq. (|^)]. In our analysis we decompose the components of the field F^i^{x) in spherical 
harmonics. Since this decomposition is nonlocal (it involves an integration over the 2-sphere r,t = const), it will 
generally depend on the extension of k'^^'^^ ^ which we now have to specify. 

A natural extension of k'^'^'^^ was prescribed by MST (also MNS (D) by setting in the right-hand side of Eq. (||) 
Ui{x) and g°'^ — > g°'^{x)\ namely, by PP the four velocity it" from zi to x, and assigning to the metric function 
its actual value at x. However, for our analysis, we found it useful to apply a different extension: one in which all 
(contravariant) tensorial components of k"^'^^{x) are assigned fixed values — the same values they have at a; zq; 

k°'<^'^\x) = kf'^^. (15) 

Note that this definition is coordinate-dependent; here we refer to (contravariant components in) the Schwarzschild 
coordinates. Throughout the rest of this manuscript, k"^^^{x) will denote specifically the extension defined in Eq. 
(|l5|), to which we shall refer as the "fixed components" extension. This extension turns out to be most convenient for 
the numerical determination of the modes of the full force (recall that the same choice of extension must be made for 
both the direct and the full forces!). 



C. Constructing the direct force 

To analyze the direct force, F^[^{x) = ijk"^^^h'^".g, we first use Eq. ( |l^ ) to obtain 

hf^.s = -2^le-^S,supu^ + 4iie-\u0,suy + upu^-s) - t-^^S^sPf^ /2 + e'^Pf^s^ (16) 

where -Pg^l^^ = P^-^-s is a regular function, of order 5x (and higher orders). Since iia is PP (from zg to a;), its covariant 
derivatives are proportional to 5x. Therefore, the second term in the above expression may be absorbed in the fourth 
term: this merely amounts to modifying the explicit form of P'p^^- Considering next the third, oc e^^ term in Eq. (|l^), 
and recalling S,s oc 0{5x), we write this term in the form e"'^^'^!^^, where Pjj^^ is a regular quantity of 0{5x^) (and 
higher orders). Absorbing then the term t'~^P'p'^j — e^'^ {^'^P'p^sj the term e"^^'^!^''^ (which amounts to re-defining 
Pf^s)^ finally write /ij^'^.^ as 

Consequently, the direct force takes the form 

Pdir(^) = {~\K'^'e-^S.s + 6-^P("3)) , (18) 

where 

j^cS ^ 4fc"/37<5y^^^^ (19) 

and nP^^^ = k'^f^^^pf^g (the factor 4 is introduced for later convenience) . Note that the quantity P^-^ , which is regular 
B,t X = zq, generally contains also terms of order Sx"^ and higher. However, the contribution from such higher-order 
terms to P^;^ vanishes at 5x = 0, and so these terms may be ignored in our analysis. We shall indeed drop these 
higher-order terms, and take P^-^ to be a polynomial in 6x of homogeneous order 5x^ . 

The coefficients of the tensor K"^ are not constant, as the field up is a PP field and not a field of constant 
components. It will prove convenient to expand K"^ in 5x, and express it as 

where Kq^[= K°'^{x Zq)] is a field of constant components, K"^ is proportional to Sx, and so on. (Note that the 
terms ^">o — unlike Kq^ — depend on the extension.) Considering now the first term in the expression for the direct 
force, Eq. (|l^), and recalling 5* ,5 cx Sx, we observe that the contribution from the term K^^ and higher order terms 
of K""^ to the direct force vanishes a.t x ^ zq. We hence drop these terms. In addition, we observe that the term 
X^^e~^S^s may be absorbed in the term e~^P^s^ of the direct force, which merely amounts to re-defining P(ly Thus, 
the direct force takes the form 
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^d"r = (^\Kfe~'S,s + lKfe-'S,s + e-'P^,^ . (21) 

As in Paper I, we now expand S in powers of 6x, in the form 

S^So + Si + S2 + ---, (22) 

where 5*0 is the leading order (oc Sx'^) term of S, Si is the correction term of homogeneous order Sx^, and so on. In 
this work we will need only the explicit form of Sq: 

So^ {g^u + u^u,)Sx''Sxr (23) 



The factor e appearing in the last expression for the direct force, Eq. (21), is then expanded as 

^-3 _ Q-3/2 _ <;.-3/2 3 5/2 / 15 7/2 2 3 5/2 



-3 3 5 „ ,-7/1^02 3 2 



7:^o'Si+e,' l-St~-et,S, +■■■, (24) 



where eg = Sq^^ . In this expansion, the first term scales as Sx~^, the second as Sx^'^, and so on. The terms included 
in the three dots scale as 6x'^ or higher powers of Sx. 

Next we expand the direct force in powers of Sx, using the above expansions of S and e^^. Based on Eq. (|p), this 
expansion takes the form 

F^-i, = /i^ (eo-3P("i) + £0 '^P("4) + 60 , (25) 

in which P^^^ denote polynomials of homogeneous order n in Sx, and where we have omitted higher-order terms 
that vanish at x zq. Notice that the term e^^P^gj of Eq. ( pT| ) has been absorbed here in the term Eq ^P^"^ (with 
higher-order corrections that vanish at a; ^ zq and are thus omitted). Also absorbed in e'^'^PIj-^ are other terms like 
Kg'^eQ^S2,5, Kg:^eQ^SiSi^5, K'^^e^^Si^s, etc. The functions P^^^ and P^l^ are given explicitly by 

=-^^o'5o,5, (26a) 
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-P(4) = ~2^o SqSi^s + -^Kq SiSo,5 — -^K" SqSq^s (26b) 

(the explicit form of P^y-j will not be needed). Note that the leading-order term of the direct force, fi^eQ^PI^^y emerges 
exclusively from the leading-order term oc K;^^ in Eq. (|l]), wher eas the next-order term, /i^e^ ^-P(4)i is composed of 
contributions coming from both terms oc Kq^ and oc iff*. The oc Kq^ contributions (and thus the entire leading-order 
term) are all analogous to ones that occur in the scalar model [see Eq. (23) of Paper I], whereas the oc iff* contribution 
has no counterpart in the scalar case considered therein. We also point out that, since if^* does not depend on the 
extension of k°'P'^^ (unlike if"*Q), one finds that the leading-order term of P^jj. is extension- independent, whereas the 
explicit form of the higher-oder terms does depend, in general, on the choice of extension. 



III. DERIVATION OF THE REGULARIZATION PARAMETERS: GRAVITATIONAL CASE 



In principle, the derivation of the RP will now involve expanding the direct force components P^jj. in scalar spherical 
harmonics, and then taking the limit a: — > zq [just as in the scalar case analysis — cf. Eq. (27) in Paper I]. This will 
yield the /-mode contribution to the direct force, P^jJ., from which one may deduce the values of all RP. However, 
at this point we may exploit the remarkable analogy between the expression derived here for the gravitation direct 
force, Eq. ([25|), and the corresponding expression obtained in the scalar toy model [see Eq. (22) of Paper I]: these 
expressions differ only in the explicit form of the three coefficients Pj." ^ . Conveniently, this analogy will now allow 
us to base most of our analysis on results already obtained in Paper I. 
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We begin by recalling the expression obtained for the direct force in the scalar case [see Eq. (27) in Paper I]0: 

j-i(dir,sca) _ 3 p(l,sca) 5 p(4,sca) 7 p(7,sca) /„„n 

where pj^''^'^^^ = —ISq^s, P^'*''''^^'' — —^SqSi^s + jSiSo,s, and pj^'^'^^^ is a polynomial in 6x, of homogeneous order 
whose explicit value will not be needed here. We use the label "sea" to distinguish quantities associated with 
the scalar case, from their gravitational-case counterparts. 

Comparing Eqs. ( p5| ) and (p7|), taking into account also the explicit form of the coefficients pj^''*'*''^'') and Pg^''^\ we 
now express the gravitational direct force as a sum of three terms, in the form 

= + F2"+F3"), (28) 

where 

rpa l^aS jTidir.sca 

^ e,' - K^'P^'^^'^^^) . (29) 
We proceed by considering separately the contributions to the RP from each of the three terms 3- 

A. Contribution to the RP from the term F" 

Consider first the term of the gravitational direct force. This term is just the scalar direct force, contracted 
with Kq^ — an array of constant coefficients [recall A'q = K^^ldx — > 0), where K"^ is the tensor defined in Eq. (p^j. 
Since the constant array Kq^ does not interfere with the multipole decomposition, one may immediately conclude 
that the contribution from the term to any of the RP, in the gravitational case, would be precisely the same as in 
the scalar case — multiplied by Kq^ . Denoting by Rf {i = 1, 2, 3) the contribution of the term F" to any of the RP, 
we thus simply have 



R'^ ^ K^^ Rf''''\ (30) 



where the scalar-case values R^f'^'^^ are those given explicitly in Paper I. (We have made here the obvious replacement 
q /i.) In particular, since '^''^ = of'"'''^ = 0, we find = P>f = 0. 



We next consider the term F2 in Eq. (Eq). This term has the from oc Cq '^P(\)- As shown in Paper I, in evaluating 



the contribution of this kind of terms to the Z-mode of the direct force at zq, one is allowed to take their limit 6t, (5r — > 
before applying the multipole decomposition. This is true regardless of the explicit form of the polynomial ^'(4). We 
hence proceed by considering F2{St = 6r = 0); we show that this quantity actually vanishes, even before applying 
the multipole decomposition. 

Examine the form of F2, as defined in Eq. ([29|): The quantity K"^, recall, is the first-order variation of the tensor 
j^aS ^ ik^'^^^UjSUj with respect to 6x. Recalling that k°'l^'>^ is a tensor of constant components, we have K"^ — 
ik"^'^^ {5uj3U^ + u/^Su^) , where Sup is the first order variation in the PP four-velocity up{x) (namely. Sup = TppUxSxP , 
with being the connection coefficients at zq). Thus, 

F^ = -2e^^k°'^'<^{Suf3U^ + upSu-y)So,s- (31) 



Consider now the explicit form of k"^^'^^, given in Eq. (^). Three of the five terms of k°'^'^^ are proportional to 
These three terms will contribute nothing to as 
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^For later convenience, we use here a re-definition of ^j<*"^'='=^)^ with the factor omitted {q is the scalar charge) 
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u^Sq^s = 2u*(.g^5 + u^us)Sx'^ = 0. 



Consider next the first and fifth terms of k"'^''^ , proportional to u'^u'' and g^^ , respectively. Both terms yield 
contributions to F2 which are proportional to 6upu^ . This quantity, in fact, vanishes for our orbital setup: To see 
that, first recall = u^up = —1, as the length of the four velocity is preserved when PP. Then, observe that the 

linear variation of this equality with respect to 5x yields 

= 5{u^up) = Su^up + ii^Sup 

= 2Su[3u'^ + Sg"^ {x)uaUf3 (32) 

(to linear order in Sx), where Sg°'^{x) = g°'^ ^{zi^)5x'^ is the linear variation in g°'^{x). Since in our setup the 
trajectory is equatorial, and since 5"^ g — g"^ ^ = at the equatorial plain, the linear variation Sg"^ vanishes (recall 
that in considering F2 we reduce dx to just 69, S(p). Consequently, we obtain from Eq. ( ^ i5m^m'^=0. 

In conclusion of the above discussion, we find F2 = (in the limit St — 5r — 0). Hence, obviously, this term yields 
no contribution to any of the RP: 

= 0. (33) 

Note that this result may no longer be valid when using fc-extensions other than the "fixed components" extension 
employed here: Usually, there will arise additional terms in Eq. (^), corresponding to first-order variations of the 
tensor k°'^^^{x). Also, notice that the result ( ^3|) will generally not hold when considering non-equatorial orbits, as 
the variation Sg""^ in Eq. ( ^ ) will generally fail to vanish. In both cases (namely, a different fc-extension, and/or a 
non-equatorial orbit), our calculation would lead, in general, to a nonvanishing contribution R^- 



C. Contribution to the RP from the term 

We finally turn to the term F^ in Eq. (|2^). Recalling that K^^ is just an array of constants, we may write this 
term as 

F3" = e-^P("7), (34) 

where P^j-^ = P^j^ — p^'^'^'^^'> jg once again a polynomial in Sx, of order Sx'^ . 

The contribution of the term F^ to the Z-mode direct force is obtained by carrying out the (Legendre) integration 
over a 2-sphere r — t — const, and then taking the limits St, Sr — » 0. As shown in Paper I, in evaluating the 
contribution of a term of the form (^J) (regardless of the explicit form of ^^(7)) one may interchange the integration 
and the limits, and set St = Sr — Q before integrating over the 2-sphere — just as with the term F2 considered above. 
To carry out the Legendre integration, it proves especially convenient — aSrin Paper I — to use a new set of spherical 
coordinates {9' , ip'), in which the particle is located at the polar axis, 9' — OB The contribution from F^ to the l-mode 
direct force can then be expressed as 

/ ^^^-f'(7)^i(cos^') d(cos 9')d^', (35) 

where Pi is the Legendre polynomial, and eo7-P(7) reductions of eoj-P(7)j respectively, to r = rp and t = t^. 

(Note that, conveniently, in the 9' , ip' system the contribution to any Z-mode at x —>■ zq comes only from the axially- 
symmetric, m = mode.) From Eq. (|2^), recalling cq = Sq^^ and ue = 0, we obtain, explicitly, 

io=[rl{59' + Sp') + ulS^'f\ (36) 

To implement the integral (^5|), it proves convenient, as in Paper I, to introduce Cartesian-like coordinates x,y on 
the 2-sphere, which we define here by 



^We should emphasize here that we do not regard 9' , Lp' as new spacetime coordinates — namely, all vectorial/tensorial quantities 
are still taken with respect to the original coordinates 6, tp. That is, 9' , p are merely used here as new variables for implementing 
the Legendre integral. The same holds for the coordinates x, y introduced below. 
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x = e'cofnp', y = e' sirup'. (37) 

Note X — y — at zq, and hence we have simply Sx^ — x and Sx^ — U- It is simple to show that a choice of 
transformation (6*, (/?) {9',(p') can be made, such that the coordinates x,y would relate to the original coordinates 
9, ip through 

x^5ip + 0{Sx^), y^59 + 0{5x'^). (38) 



Expressed in terms of the new coordinates, the polynomial P^"^ {69, 5ip) in Eq. ( |35| ) becomes P^j-^ {y, x) + 0{Sx°), where 
P^"-) is a polynomial of homogeneous order 7 in x, y. The contribution from the 0{6x^) corrections to the direct force 
vanishes at x — > zq: a-nd can therefore be omitted. From Eq. ( ^6| ) we also get 

eo = eoix,y) + Oi6x^), (39) 

1 /2 

where eo = [''ol^^ + 2/^) + u'^x'^] ■ Again, only the leading-order term here contributes to the direct force at a; — > zq, 
and we are allowed to drop the 0{Sx^) correction. Hence, as far as the calculation of the RP is concerned, we may 
express the contribution from the term Pf to the Z-mode direct force, Eq. (^5|), as 



/ ^o'P(7)i^.y)Piicos9')dxdy. (40) 

Note that the Jacobian of the transformation (cos0', (p') {x, y), which is actually given by 1 + 0((5x^), has been set 
here to just 1: The higher-order corrections are once again omitted, as they vanish at x —> zq- 

Examine now the integral in Eq. (^): Eq is an even function of both x and y. So is the function cos9' . However, 
each of the possible individual terms in the polynomial P^-^ (like cx xy^, or cx x'^y^, for instance) is necessarily an 
odd function of either x or y. Consequently, we observe that the entire integrand in Eq. ( |40| ) is odd in either x or y. 
Therefore, obviously, the integration over the 2-sphere would vanish. As a consequence, no contribution to the RP 



will arise from the term 



^3 = 0. (41) 



Notice that this last result is valid for any (sufficiently regular) fc-extension. A modification of the extension would 
only affect the explicit form of the polynomials P^y-j , but would not alter the odd parity structure of the integrand in 
Eq. (P. 



D. Summary: RP values in the gravitational case. 

Let us now collect the above results: We have found that neither of the terms P^ and P^ contributes to the Z-mode 
direct force. The sole contribution to the RP comes from the term P" — this contribution is given in Eq. (pO|). The 
RP in the gravitational case are therefore given by 

R"" ^ K^^ Rf''''\ (42) 

where, recall, R" stands for any of the RP, and the scalar-case values pf'^'^'^ are those given explicitly in Paper I. We 
now need only to provide the explicit form of Kq^: Recalling Kq^ = K"^{x — > zq), one easily gets from Eqs. ( |l9| ) and 

Kf^g°'^^u'^u\ (43) 

where, recall, and g""' denote the values of these quantities at zq- Note that K^^ is just the spatial projection 
operator at zq, namely, (^Kq^Vs) = for any vector Vg. 

Let us finally write Eq. ( ^^ more explicitly: First, recalling (see Paper I) that the scalar parameter A'f'^^^ has no 

component tangent to (i-e., A^^'^^^ =0), we simply obtain 

= ^fsca) (44) 

(with the obvious substitution g ^ /i). Unlike the situation with the parameter A", the quantity u^Bg^'^^^ does not 
vanish [see Eq. (85) in Paper I], and we leave the expression for B" in the form 
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= K^^Bf'"'^ (45) 

(again, with q Finally, as c'f'^'^^ — uf'^^^ — 0, we shall have, in the gravitational case alike, 

C" = £)" = 0. (46) 

IV. DERIVATION OF THE REGULARIZATION PARAMETERS: ELECTROMAGNETIC CASE 

In this section we consider the electromagnetic self force acting on an electrically charged particle: we prescribe 
the mode sum scheme in this case, and construct all required RP for an arbitrary (equatorial) geodesic orbit in 
Schwarzschild spacetime. The same analytic calculation used for deriving the gravitational-case RP will prove directly 
applicable also to the electromagnetic case, with only minor adaptations required. 

We shall consider a particle carrying an electric charge e (with |e| <C M), and assume the same orbital configuration 
as in the gravitational case (namely, the particle is taken to move along an equatorial orbit, which in the limit e — > 
becomes a geodesic). We shall also maintain here the notation for the various quantities zq, x, e, S, Ua, and so on. 
We shall denote by (pa (x) the vector potential associated with the "full" electromagnetic field induced by the particle. 
In this section we ignore the gravitational self force. 

A formal expression for the electromagnetic self force in curved spacetime was obtained long ago by DeWitt and 



Brehme 15 1 (and was reproduced recently by QW ^ using a different method). For a geodesic in vacuum spacetime. 



the electromagnetic self force is obtained from an electromagnetic "tail" force, just as in the gravitational case: 

= lim Ffjr^(-), (47) 



where hereafter we use the label "EM" to signify quantities associated with the electromagnetic case. The formal 
construction of the vector field F"^^^\x) is described in jl^,|^. Like in the gravitational case, the electromagnetic 
tail force can be written as the difference between a "full" force and a "direct" force — just as in Eq. (|). In the 
electromagnetic case, these two vector fields are given by ^ 

where (f>'^"{x) is the "direct" part of the vector potential (given explicitly below), and k'^^''{x) is a (sufficiently regular) 
extension of the tensor 

kf'^=g"^u'^-g"^u^, (49) 

defined at zq. As in the gravitational case, we shall adopt here the "fixed components" extension, defined (in 
Schwarzschild coordinates) through k°'^'^{x) = k'^^'^ . 

The mode sum prescription for the electromagnetic self force is completely analogous to the one prescribed in the 
gravitational and scalar cases: Given the (scalar harmonic) Z-modes F^-^^^'^ of the electromagnetic full force, the 
electromagnetic self force is constructed through 



OC 

pQ(EM) _ 

(=0 



hm F"T^\x) - A"(EM)^ „ 5a(EM) _ ^a(EM) 

x^za 



jjaiEM)^ (50) 



where the various electromagnetic-case RP are to be obtained, again, by analyzing the niultipole modes of the direct 
force. The rest of this section is devoted to calculating these electromagnetic RP. 

As in the gravitational case, our starting point would be the expression for the direct part of the particle's field — this 
time the direct part of the vector potential — as obtained by MNS [see Eq. (B3) of Ref. In precisely the same 
manner as in the gravitational case, this expression can be brought to the form 



t>f{x) = eS-^/^iifj + S-^'''pf{x) (51) 



[in analogy with Eq. (|14|)], where Pg is a (regular) function of 0{5x ). The derivatives of the direct vector potential 
then take the form 

cpf^ = -{e/2)S-^'^S,,up + eS-^'^up,, ~ S'^'^ S^.Pf^ 12 + S-^'^Pf}^, (52) 
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which [as in obtaining Eq. (|17|)] can be put into the form 

(^^f, = -iee-35,>^ + e-3p(3\ (53) 

with Pj^^ being a (regular) function of order 6x'^ (and higher orders) . Consequently, the direct electromagnetic force 
takes precisely the same form as in Eq. (p^, 

where this time 



F,r"i^) = [ -^K^M^-'S.s + 6-^P("3) ) , (54) 



^(EM) ^ (55) 

and ePg) = k^l^-^pf^. Again, we may drop all terms of P^^^ which are of order 5x^ and higher (as they do not 
contribute to the direct force at a; — + zo), and take P^g^ to be of homogeneous order 5x^ . 

Thanks to the complete analogy between the forms of the electromagnetic and gravitational direct forces [compare 
Eqs. (^) and (|lj)], our analysis now proceeds precisely as in the gravitational case: We expand K'^^^-^ and S in 5x^ 
as in Eqs. (|2^) and (p2|), and consequently write the direct force as a sum of three terms — in precise analogy with 
Eqs. (H) and (H): 

pa(EM) _ 2 / -3 pa(EM) 5 pa(EM) 7 pa(EM)N 

^dir — ^ y'O -^(1) '0 -^(4) ^ ^0 ^(7) ) 

^ e2 ^po(EM) ^ ^a(EM) ^ ^o(EM)^J ^ ^^g^ 

where ^"2^'^'' defined in Eq. (|29|), with the replacements if"'' — > and P^"^ ^ Pj^j^^''. The only point at 

which our current analysis differs from the gravitational case is in the explicit values taken by the various coefficients 
K^^ (and consequently in the explicit values of the terms -P("))- 

Consider first the contribution to the Z-mode direct force coming from the term i^"*^^^-': In a complete analogy 
with Eq. (pO|), we obtain 

^?(EM)=^ofEM)4'''^('?-e), (57) 

where ^"(em) (* ~ ^' stands for the contribution of the term p°''^^^^ to any of the RP, and the array of constant 
coefficients K^Ly^-^ = Kr^y^Ax — > zq) is now given by 



Noticing K^L^^-^ = -K^^ [compare Eq. (|8|) to Eq. (O)] and recalhng Eqs. (p^ and (E|), we then conclude 



^o(EM) = -(.<?"' (58) 

0(EM) 

^?(EM) = -^?(grav) (M - e) = -Pfg,,,) {^i^e), (59) 

where hereafter we use the label "grav" to signify the gravitational-case values. 

Next, consider the term i^^^^^-' ~ I^i[-em)^o,s ■ Here, the coefficient K'^^-^y^-^ (the first order correction in 

■^(EM)) is given by 

^UEM) = k^^'5up. (60) 

As in the gravitational case, it is easy to show that p!^''^^^ (evaluated at 5t = 5r = 0) actually vanishes, even before 
taking its multipole decomposition: From Eq. ( p9| ) we observe that k^-f^^ is composed of two terms, one proportional 
to and the other proportional to . The cx iT" term contributes nothing to ^^"(em), since u^Sup = (as explained 

when discussing the gravitational case). The oc term will yield a zero contribution as well, by virtue of So,su^ = 0. 

We thus find that in the electromagnetic case — ^just as in the gravitational case — the term p^^-^^'^ contributes nothing 



to the RP, namely, -R2(em) ~ ^ 
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As to the last term in the electromagnetic direct force, F^^^^^: Using the same parity considerations as in the 
gravitational case, one shows that the contribution from this term to any of the RP will vanish, namely, ^3(gM) = 0. 

This vanishing is irrespective of the explicit form of the polynomial P^"^^'^-' in Eq. (|56|). 

In conclusion, thus, we find that the sole contribution to the RP in the electromagnetic case comes from the term 
^a(EM)^ Eq. (H). We hence obtain 

^"em) = -R{gia.v)if^ ^ e). (61) 



V. SUMMARY AND CONCLUDING REMARKS 

Let us summarize our mode-sum prescription for constructing the gravitational and electromagnetic self forces. We 
start with the gravitational case: 

1. For a given trajectory, compute the tensor-harmonic modes of the metric perturbation, /i^*)' ™(r, i), by numeri- 
cally integrating the separable field equations (e.g., in the harmonic gauge [Q). 

2. Given /i^*)' ™(r, t), construct the full modes Ff^'j at the particle's location. This is done by applying the op- 
erator in Eq. (||) to /i^*'' ™(r, t), using the "fixed components" extension described above, and then expanding 
the resultant field into scalar spherical harmonics (and summing over il'm for a given I). This procedure is 
implemented in Refs. [^,0. 

3. Use Eqs. ( |6^ ) below [along with Eqs. (83) of Paper 1] to obtain the RP values corresponding to the trajectory 
under consideration. 



4. Finally, apply the mode-sum formula, Eq. ( |11| 

This prescription is now being implemented by Barack and Lousto for radial and circular orbits in Schwarzschild 
spacetime. 

The prescription for constructing the electromagnetic self force is similar: First, one has to compute the vector- 
harmonic modes of the (full-field) vector potential for the given orbital configuration. Then, one constructs the full- 
force modes Pf"|P'^^ — this construction is carried out by applying Eq. (^8|) to each of the full-field vector- harmonic 
modes, and then decomposing each of these modes in scalar spherical harmonics. Finally, one applies the mode-sum 



formula (50), with the electromagnetic RP values given in Eq. (k32) below 



The values of the RP in the gravitational and electromagnetic cases are summarized as follows: 

^(grav) ^ _^(EM) ^ ^(sca) ^ (g2a) 
^(grav) ^ _5(EM) ^ ^ ^(sca) ^ ^g2^^ 

^(grav.EM) _ jj{giav,EM) _ g (62c) 

(with the obvious replacements q ^ ^ or q ^ e), where the quantities labeled "sea" are the scalar-field parameters 
given explicitly in Eqs. (83) of paper 1. The gravitational RP were calculated previously in a different method, using 
the Z-mode Green's function expansion technique |§,D, in the special case of radial orbits. The results agree with the 



values (62) 



As we discussed above, there is a certain ambiguity in the values of the RP, which arises from the freedom in 
choosing the extension of the tensor k off the evaluation point zq (the choice of this extension affects, of course, the 
multipole decomposition of the force). However, our mode-sum scheme produces no ambiguity in the eventual value 
of the self force — one only has to make sure that the full-force modes F^^^ in Eq. (^l]) [or in Eq. (^0|)] are calculated 
using the same extension as the one used in calculating the RP. It is thus essential to recall here that the RP values 
summarized above are those referring to the "fixed components" extension of the tensors k"^^^ or k"^^ (expressed 
in Schwarzschild coordinates). This extension is most easily applicable in the numerical computation of the full- force 
modes mil. 

Based on our above analysis, we may phrase the following general statements concerning the extension-dependence 
of the RP in the electromagnetic and gravitational cases: (i) The RP A", C", and are insensitive to the extension 
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of k°'f^'^^ (provided it is regular enough), (ii) The value of does depend, in general, on the choice of extension; 
however, all sufficiently regular extensions which differ from the "fixed com ponents" extension k°'Pi^ by an amount 
of only 0{5x^) will admit the same value of i?" — the one given in Eq. (62b). It is interesting to refer here to MNS's 



analysis [||, in which a different extension has been employed: MNS extended the tensor k"^'^^ by PP the four velocity 
from zq to X and just assigning to g"^ the actual value it has at x. Interestingly, within this extension [differing from 
the "fixed components" extension already at 0{Sx)], all RP attain precisely the same values as in the scalar case |^ 
(except that in the electromagnetic case all RP are to be multiplied by —1). 

Finally, it is important to remind that the gravitational self force is a gauge-dependent notion, as discussed in Ref. 
1^. The prescription described in this manuscript applies to the self force associated with the harmonic gauge (in 
which the original MST/QW scheme has been formulated). It also applies, with the same RP values, to any other 
gauge related to the harmonic gauge through a regular gauge transformation However, for other, non regular 
gauges, the mode-sum scheme is not guaranteed to be valid in its present form. A method for overcoming this gauge 
problem has been sketched in and is currently being implemented for circular orbits in Schwarzschild p^ . A 
different strategy (applicable in the Schwarzschild case) would be to calculate the self force directly in the harmonic 
gauge 
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APPENDIX A: DIRECT METRIC PERTURBATION 

In this appendix we obtain Eq. (^3|) for the direct metric perturbation, by processing the expression given by MST 
in Ref. §. 

By considering the Hadamard expansion of the (full) metric perturbation, MST obtained the following expression 
for the (retarded) trace-reversed perturbation [sec Eq. (2.27) of Ref. §]: 



hf3j{x) = 2/i 



2l3{x,zo) 



tail term + O(e^). (Al) 



Here we use the notation of our Sec. ( see Fig. |l|), namely: x is a point in the neighborhood of the force evaluation 
point Zq] e = S^/^ is the length of the short geodesic section connecting x to the worldline and normal to it; zi 
denotes the intersection of this geodesic with the worldline; and (or ul^) denote the four- velocities at zq and 
zi, respectively; and u°' and iif (or u\) arc their PP to x. In addition, parentesized indices denote symmetrization, 
and Rj^/^^s represents the Riemann tensor PP from zi to x with respect to any of its indices carrying the hat sign. 
The "tail term" represents a non-local contribution to the full perturbation, with its form given explicitly in The 
function (3 (denoted k"^ in [Q) is a regular function satisfying (3—1 + 0{5x^) [see Eq. (A14) therein]. Note that the 
correction term proportional to the 4-acceleration in Eq. (2.27) of Ref. [|l|] can be omitted, since, for geodesic orbits, 
it contributes to the self force only at order higher than 0(/i^). For the same reason, we omit here the 0(r~^e) term 
indicated therein. Finally, notice the notational change a 5/2. 

The direct part of the metric perturbation is now taken as the difference between the full perturbation given in Eq. 



(Al), and the tail term. The terms included in 0(e ) do not contribute to the direct perturbation at the limit x ^ zq 
(e ^ 0); nor do they contribute to the direct force, whose construction involves only first-order derivatives of haf){x). 
We thus re-define the direct perturbation by ignoring these O(e^) terms, which leaves us with only the three terms in 
the squared brackets, scaling as e~^, e^, and e^, respectively. 

Consider now the second and third terms in the squared brackets: First, note that the (coordinate components of 
the) two vectors u ? a nd differ only at 0{5x). Hence [recalling e,S''^ oc 0{5x)], this difference contributes only 
to 0{5x'^) in Eq. (^). We may thus ignore this correction, and just replace with u" in the second and third 
terms in the squared brackets. Likewise, we replace u" with in the second term. Similarly, we may ignore the 
0(5x) difference between R]^jj^g and Rap-yS (the latter denoting the coordinate-value of the Riemann tensor at zq) as 



it contributes only to 0{6x^) in Eq. (Al). The direct metric perturbation thus takes the form 

hf:^{x) = 2^. 



(A2) 
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Examine now more closely the second term in Eq. (A2): Since S'"'^ cx 0{6x), the only contribution to the direct 
force which does not vanish at x ^ zq arises from differentiating S'''^. RecaUing Eq. (^), we have 



u us) 



0{Sx) 



Note that S'^^u^ 



(at a; ^ zq). Note also that the second term in Eq. (A2) (unlike the other two terms) is perfectly 
regular at x = zq- This allows us to evaluate its contribution directly at zq, which we do by just "removing" the hat 

Recalling Eq. (0), the contribution from this term to the direct force at zq then reads 



symbols from R^xpii ^"^^ 



(evaluates at a; = zq)- Examining the form of the tensor k"^'*^ , given in Eq. (H), we observe that three of its five terms 



are proportional to u , and thus vanish when contracted with S' 



The first term of k"^'^^ is proportional to u^u^ , 
and thus yields a vanishing contribution when contracted with either u'^RaXpp or u^Ra^xpj. Likewise, the last term of 
l^ap-)S^ proportional to g''^, is found to vanish when contracted with either uPRaXp/sUj or u^RaXp^up. We conclude 
that the contribution of this regular term to the self force vanishes — even before taking the harmonic decomposition. 

Finally, consider the third term in Eq. ( A2). Noticing that this term has the form ex (a regular function of x), we 
may write it as 



(A3) 



where Pp-* is a certain regular function, of order 5x^ . (This form will be convenient for our analysis in Sec. II.) We 
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further notice that the terms of 0{5x'^) included in the function [i contribute to hp^ an amount of precisely the form 
(A3). We may thus absorb this contribution in the contribution (A3) coming from the third term, and replace the 
function (3 with just 1. The explicit form of the regular function P^^ {x) will not be needed in our analysis. 

In conclusion, we find that the direct part of the metric perturbation is effectively given by Eq. (^3|) . This expression 
is used as a starting point for the analysis in this paper. 



[1] Y. Mino, M. Sasaki, and T. Tanaka, Phys. Rev. D 55, 3457 (1997). 
[2] T. C. Quinn and R. M. Wald, Phys. Rev. D 56, 3381 (1997). 

[3] L. Barack and A. Ori, Phys. Rev. D 61, 061502(R) (2000); L. Barack, Phys. Rev. D 62, 084027 (2000). 
[4] L. Barack, Phys. Rev. D 64, 084021 (2001). 

[5] L. Barack and C. O. Lousto, Phys. Rev. D ??, ??????(R) (2002), |gr-qc/020504^ . 

[6] L. Barack, Y. Mino, H. Nakano, A. Ori, and M. Sasaki, Phys. Rev. Le tt. 88, 091101 (2002). 

[7] L. Barack and A. Ori, Phys. Rev. D ??, ?????? ( 2002), ^r-qc/0204093| . To be referred to in this manuscript as "Paper I". 
[8] Y. Mino, H. Nakano, and M. Sasaki (submitted), |gr-qc /Oil 1074 
[9] L. Barack and A. Ori, Phys. Rev. D 64, 124003 (2001). 
[10] An alternative, elegant formulation of the self force was recently introduced in S. Detweiler and B. F. Whiting (submitted), 
|gr-qc/ 0202084 

[11] Y. Mino and H. Nakano, Prog. Theore. Phys. 100, 507 (1998). 

[12] A closely related calculation scheme, also employing a mode sum decomposition (but relying on a different regularization 
scheme, based on the zeta-function regularization technique), was introduced in C. O. Lousto, Phys. Rev. Lett. 84, 5251 
(2000); 

[13] A full implementation of the mode-sum scheme for the scalar self force is demonstrated in L. M. Burko, Phys. Rev. Lett. 84, 
4529 (2000); L. Barack and L. M. Burko, Phys. Rev. D62, 084040 (2000); L. M. Burko, Y.-T. Liu, and Y. Soen, ibid. 63, 
024015 (2001); and L. M. Burko and Y.-T. Liu, ihid. 64, 024006 (2001). The mode-sum scheme has also been implemented 
in the gravitational case, for radial Q and circular |l^ trajectories around a Schwarzschild black hole. 

[14] L. Barack and C. O. Lousto (work in progress). 

[15] B. S. DeWitt and R. W. Brehme, Ann. of Phys. 9, 220 (1960). 

[16] N. Sago, H. Nakano, and M. Sasaki (submitted), ^r-qc/020806c| . 



14 



